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1 Introduction 

1.1 Background about Wavelets 

Wavelets are functions used to approximate other functions or data. They are particulary effective 
in approximating functions with discontinuites or sharp changes. Wavelets further preserve both 
the global picture and local detail of a function. Applications of wavelets include compressing 
images, such as X-rays and the FBI's fingerprint collection, and denoising data, found in many 
areas, including geology, meteorology, astronomy, acoustics, and econometrics. 

1.2 Discussion of 

D4 is the univariate Daubechies scaling function with four refinement coeffiecients. D4 may be used 
to exactly reproduce univariate functions which are constant or linear. 

2 Overview of Bivariate Daubechies Scaling Function 
2.1 Conditions on (f){x,y) 

The idea of this project is to create scaling functions y) of two variables which satisfy properties 
of the type satisfied by The dilation equation for a bivariate Daubechies scaling function [D] 
will be: 

4>{x,y) = '^'^Cij(l){2x -i,2y - j). (1) 

* j 

where the sums are over all integers. This scaling function will satisfy several conditions. One 
condition to satisfy is compact support. The support of the bivariate scaling function is from to 3 
in the x direction and from to 3 in the y direction. The function (p is zero outside of this support, 
and is also zero on the boundary. There is also an averaging condition that comes from the idea 
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that the average value of the scaling function is not 0. Related to this condition is the assumption 
that the inner product 

/oo roo 
/ (f>{x,yfdxdy = 1 (2) 
-oo J — oo 

A bivariate scaling function must be orthogonal to its three types of integer translates: those in 
the X direction only, those in the y direction only, and those in both the x and y directions. This is 
called the orthogonality condition. (The orthogonality is important in order that the function and 
its translates form a basis.) Also, the scaling function must be designed to exactly reproduce linear 
functions of the form f{x, y) = kx + ly + ni, where k, I, and m are constants. This last condition is 
called the regularity condition, and it also guarantees that (f) is continuous. (However, Han points 
out that any (j) we create will not be differentiable everywhere, a property shared by -D4 [H].) The 
consequences of all four conditions will be explored in the next several sections. 



2.2 A Formula for the Refinement Coefficients 



Due to the orthogonality condition, the set of functions 

{0(2x-i,2y-j)|i,j GZ} (3) 

is a basis of a vector space of which (pix, y) is a member. This observations allows us to prove this 
Lemma: 

Lemma 1: For all integers i and j, Cij = 4 < y), (j){2x — i,2y — j) >. 
Proof of Lemma 1: According to the Orthogonal Decomposition Theorem, 

< (/)(2x -i,2y- j),(t){2x -i,2y- j)> 

and that 

^ < 4>{x,y),(pi2x - i,2y - j) > 

< (t>{2x -i,2y- j)cP{2x - i,2y - j) > ^ > 

The denominator in this formula turns out to equal \. To see this, observe that the definition 
of inner product, along with (|2|) yields 

00 fOO 



< (p{2x - i,2y - j)(p{2x -i,2y - j) > = / (j){2x - i,2y - j)(j){2x - i,2y - j)dxdy 

J — CXD J ~00 

poo poo 

= 4/ / Hu-i,v-j)(piu-i,v-j)dudv 
1 



—00 J —00 



^ < <P{x,y),<t>{x,y) > . 

Therefore, 

< (/.(2X -i,2y- j)ct>{2x -i,2y- j) >= ^, (6) 
and we have proven that q = 4 < y), (j){2x — i,2y — j) >. 
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2.3 Consequence of Compact Support: Only 16 Non-Zero Refinement Coeffi- 
cients 

Next, assuming that the compact support for a bivariate Daubechies scahng function is < x < 3 
and < y < 3, we prove this lemma: 

Lemma 2: For i / 0, 1, 2, 3 or j / 0, 1, 2, 3, Cij = 0. 

Proof of Lemma 2: From the definition of inner product, we know that 

/oo poo 
/ <Pix,y)(l){2x - i,2y - j)dxdy. 
-oo J —oo 

When i < —2 or i > 6, and j < —2 or j > 6, then (j){2x — i,2y — j) will be non-zero where 
(j){x, y) = 0, and vice versa. Therefore, we can definitely conclude that Cij = unless —2 < i < 5 
and -2<j< 5. 

It remains to show that Cij = for i = —2, —1,4, or 5, and j = —2, —1, 4, or 5. 

We can use our knowledge that some refinement coefficients are zero in order to prove that other 
ones are zero. The key to showing that the other coefficients are zero is the following observation: 
for any integers m and n, 

TYl Tl 

<A(y , 2 ^ Cm-I,n-1</'(1, 1) + Cm_l,„_2</'(1, 2) + Cm_2,n-l</'(2, 1) + Cm-2,n-2(p{2, 2) (7) 

due to the compact support requirement. For instance, 

-^) = C-2,-20(l, 1) + C_2,-30(l, 2) + C_2,-30(2, 1) + C_3,_30(2, 2). (8) 

Since any refinement coefficient which has -3 as one of its indices is already known to be zero, and 
due to the compact support condition, (^) becomes 

= C_2,-2</'(l,l). (9) 

Already we can conclude that c_2,-2 = or 1) = 0. 

From an assumption about one of 0(1, 1), 0(1, 2), 0(2, 1) or 0(2, 2) being non-zero, we can prove 
that more than half of the remaining refinement coefficients that we are trying to eliminate are, in 
fact, zero. We will demonstate how this works if we assume 0(1, 1) ^ 0. 

If 0(1, 1) 7^ 0, then c_2,-2 = 0. Then using m = —1 and n = 0, (^) becomes 

0(-^, 0) = C_2,-10(l, 1) + C„2,-20(l, 2) + C_3,-10(2, 1) + C_3,-20(2, 2). (10) 

So, by previous results, = c_2,-i0(l, 1), and therefore c_2,-i = 0. In a similar manner, it can be 
shown that c„2,j = for all remaining j, and nearly identical arguments can be made to show that 
Ci-2 = for all remaining i. 

Refinement coefficients of the form Ci-i and c^ij can be disposed of in the same way. For 
example, using m = and n = — 1, (0) becomes 

Tl 

0(0, --) = C_i,_20(l, 1) + C_l,_30(l, 2) + C„2,-20(2, 1) + C_2,-30(2, 2) (11) 

Again, by previous results, this equation reduces to = c_i,_20(l, 1), so c_i,_2 = 0. Other 
coefficients are then dealt with in the same way. Consequently, we can eliminate more than half of 
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the refinement coefficients that we are trying to show are zero by an assumption that (p{l, 1) 7^ 0. 
The same can be said of the other three values of (p. 

In the case where both ^(1, 1) and 0(2, 2) are non-zero, then all of the unwanted coefficients are 
clearly eliminated. The same is true if both (/>(!, 2) and (j)(2, 1) are non-zero. 

Without loss of generality, if (p{l, 1) and (p{2, 1) are both zero, and the other two are non-zero, 
then a few more steps are required to reach our result. In this case, there are eight refinement 
coefficients that are not determined to be zero by our previous arguments: Cij where i can be 0, 
1, 2, or 3, and j is —2 or —1. However, the same type of argument can be used to eliminate these 
one at a time. For instance, using m = 1 and n = 0, (0) becomes 



The first and third terms are zero by our assumption in this case. The coefficient C-1-2 is zero 
follows from the assumption that (/>(!, 2) is non-zero, using our arguments above. Therefore, we 
have that = co-2(p{l, 2), and we can conclude that co,-2 = 0. All eight can be eliminated in this 
way. 

Finally, if three of the four function values are zero, we can use the fourth value to eliminate 
more than half of the coefficients and then arguments such as in the previous paragraph to show 
that the rest are zero. We omit the details. 

If all four of the function values are zero, then (p is the zero function. 

In summary, for all cases, all of the refinement coefficients are zero unless < i < 3 and 
< J < 3. 

3 Deriving the 14 equations 

3.1 The Averaging Equation 

In this section we prove that 






(13) 



Here we will use the averaging assumption, which states that 




(14) 



Combining this with (|^) and (P), we can calculate the following: 



/ / (t){x,y)dxdy 




(^ ^ Ci,j(t)(2x -i,2y- j))dxdy 
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Note that the last step is vahd because the sums are finite. As a result we have the following 
equation: 

/ (l){x,y)dxdy = j^^^Cij / (p{x,y)dxdy 

oo J —etc ^ • • J ~oo J ~oo 

^ J 

Due to our averaging condition, we can divide both sides by the double integral, leaving us with 
1 = i Ei Ej or 4 = ^. Cij. 

3.2 The Orthogonality Equations I 

By the orthogonality condition, we assume that translates in the x or y direction (or both) are 
orthogonal to the scaling function. Using b and d to translate the variables, the dilation equation 
dl]) becomes: 

^(x - 6, y - d) = ^ ^ Cij(/)(2(x -b) -i, 2{y - d) - j). 

i j 

Letting m = 2b + i and n = 2d + j leads to 

(/){x - b,y - d) Cm~2b,n~2d<P{'^x -m,2y -n). 

m n 

We are interested in the inner product of (j){x, y) and (f){x — b,y — d). That is, 
< Hx,y),Hx -b,y -d) > 

= < X] Cij(p{2x -i,2y - j), ^ ^ Cm-2b,n-2d4>i'^x - m,2y - n) > 

i j m n 

= X] X] CijCm-2b,n-2d < 4'{2x - i,2y - j), (t){2x - m,2y - n) > 

i j m n 

" ^ XI XI X] XI CijCm_26,n-2d <4>{x-i,y - j), (p{x -m,y-n) > 

i j m n 

with the last line being true due to (^). 

When i ^ m or j ^ n, < (f){x — i,y — j), 4>{x — m,y — n) >= 0, by the orthogonality condition. 
So, in the previous paragraph, Ylm Yin ^^^y counts when i = m and j = n, and it is actually the 
case that < (p{x — i,y — j), 4>{x — m,y — n) >= 1. So we have that 

< 4'{X, y), </>(x - fc, 2/ - C^) >= X] X] Ci,iCi-2fe,j-2d- (15) 

i j 

If b and d are not both zero, then 

= X^ X^ CijCm-2b,n-2d- (16) 
i 3 

Specific values of b and d will lead to specific orthogonality equations. These equations take into 
account Lemma 2, which showed that only 16 refinement coefficients are non-zero. If b and d are 
both 1, then (|l6|) becomes 

C3,3Cl,l + C3,2Cl,o + C2,3Co,l + C2,2Co,0 = (17) 
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Similarly, if 6 7^ and d = 0, or if 6 = and d 7^ 0, respectively, ( |l^ ) reduces to 

C3,3Cl,3 + C3,2Cl,2 + C2,3Co,3 + C2,2Co,2 + C3,lCi,i + C3,oCl,o + C2,lCo,l + C2,oCo,0 = (18) 

and 

C3,3C3,1 + C3,2C3,o + C2,3C2,1 + C2,2C2,0 + Cl,3Cl,l + Ci,2Cl,0 + Co,3Co,l + Co,2Co,0 = (19) 

No other orthogonality equations can be derived from ([l6[). 



3.3 The Orthogonality Equations II 

There is one more equation to be derived based on the orthogonality condition. Applying the 
dilation equation (|l|) to (§) gives us: 

< 4>{x,y),4>{x,y) > = <^^Cij(t){2x -i,2y - j),^^Ckj(l){2x - k,2y -I) > 

i j k I 

= ^^^^<^ij(^k,l < - i,2y - j),4'{2x - k,2y - I) > 

i j k I 

= < "t^C^^ -h'^y- j), 0(2a: - i,2y - j) > 

i j 



<(t>{x-i,y- j), (t>{x -i,y-j)> 



4 

« J 



with the last two lines being true by the same type of argument that led to (16). Since 

< (t^ix - i,y - j)A{x - i,y - j) >=< (t){x,y),(t){x,y) >, 

we have shown that 



Due to (|2|), < (/)(x, y), (/)(x, y) > can be divided out of both sides of (pO|), leaving ^ cf j = 1, 



or 



EE^L=4- (21) 

i j 

Due to Lemma 2, the sum in this last equation is made up of 16 terms. 
3.4 The Regularity Equations 

In [J] Jia states that, for a positive integer k, a sequence a on satisfies the sum rules of order k 
if 

a(2/5 + e)p{2p + e) = ^ a{2(3)p{2p) Ve G G n^.i, 

where n^.i is the set of polynomials with total degree less than k and 0, is the set of vertices of the 
unit cube [0, 1]'^. Jia proves that satisfying the sum rules is equivalent to satisfying the regularity 



condition. We can apply the result of Jia to the case where s = 2 (the bivariate case) and k = 2 
(linear functions are reproduced). Then nine regularity equations follow from the sum rules. 

When p = l and e = (0, 1), C2i,2j+\ = C2i,2j, or 

C0,1 + Co,3 + C2,l + C2,3 = Co,0 + Co,2 + C2,0 + C2,2 (22) 

When p = l and e = (1, 1), X) C2i+i,2j+i = Yl C2i,2j, or 

Cl,l + Ci,3 + C3,i + C3,3 = Co,0 + C0,2 + C2,0 + C2,2 (23) 

When p=l and e = (1, 0), X) C2i+i,2j = Yl (^2i,2j, or 

Cl,0 + Cl,2 + C3,o + C3,2 = Co,0 + Co,2 + C2,0 + C2,2 (24) 

When p = X and e = (0, 1), Y 2«C2i,2j+i = 2ic2i,2j, or 

2C2,0 + 2C2,2 = 2C2,1 + 2C2,3 (25) 

When p = X and e = (1, 1), Y'^^ + C2j+i,2j+i = Y 2iC2j,2j, or 

Cl,l + Cl,3 + 3C3,l + 3C3,3 = 2C2,1 + 2C2,3 (26) 

When p = X and e = (1, 0), X] 2i + C2i+i,2j = Y '^ic2i,2jV., or 

Cl,0 + Cl,2 + 3C3,o + 3C3,2 = 2C2,1 + 2C2,3 (27) 

When p = y and e = (0, 1), Y 2j + C2j,2j+i = 2jc2i,2j, or 

C0,1 + C2,l + 3co,3 + 3C2,3 = 2co,2 + 2C2,2 (28) 

When p = y and e = (1, 1), Y + C2i+i,2j+i = Y '^jc2i,2j, or 

Cl,l + C3,i + 3ci,3 + 3C3,3 = 2co,2 + 2C2,2 (29) 

When p = y and e = (1, 0), Y ^3C2i+i,2j = Y '^jc2i,2j, or 

2C1,2 + 2C3,2 = 2co,2 + 2C2,2 (30) 

Note that all of these equations are linear in the refinement coefficients. 

4 Solution of the 14 equations 

4.1 Method of Solution of the Fourteen Equations 

We now have fourteen equations with sixteen variables. To review, here are the fourteen equations: 
Averaging equation: 



Orthogonality equations: 

C3,3Cl,l + C3,2Cl,0 + C2,3C0,1 + C2,2Co,0 = 
C3,3Cl,3 + C3,2Cl,2 + C3,lCl,l + C3,oCl,3 + C2,3C0,3 + C2,2Co,2 + C2,lCo,l + C2,oCo,0 = 
C0,3C0,1 + Cl,3Ci,i + C2,3C2,1 + 03,303,1 + Co,2Co,0 + Cl,2Cl,0 + C2,2C2,0 + C3,2C3,o = 

Regularity Equations: 

C0,1 + Co,3 + C2,l + C2,3 - CQfl " Co,2 - C2,0 - C2,2 = 

Cl,0 + C3,0 + Cl,2 + C3,2 " Cq.O " Co,2 " C2,0 " C2,2 = 
Cl,l + Ci,3 + C3J + C3,3 - Co,o - Co,2 - C2,0 - C2,2 = 
2C2,0 + 2C2,2 - 2C2,1 - 2C2,3 = 
2ci,2 + 2C3,2 - 2co,2 - 2C2,2 = 
Cl,0 + Cl,2 + 3C3,0 + 3C3,2 - 2c2,0 - 2C2,2 = 
Cl,l + Ci,3 + 3C3,1 + 3C3,3 - 2C2,0 - 2C2,2 = 
C0,1 + C2,l + 3co,3 + 3C2,3 - 2co,2 - 2C2,2 = 
Cl,l + C3,l + 3ci,3 + 3C3,3 - 2co,2 - 2C2,2 = 

Due to the complexity of the problem, Maple is unable to solve the system of fourteen equations 
directly. As a result, we had to take another approach. We began by observing that the regularity 
equations and the averaging equation form a non-honiogcncous system of linear equations. Solving 
this system resulted in solutions for the ten following variables ci,2, C2,i, C2,o, co,2, co,i, ci,o, ci,i, 03,0, co,3, 
and Co,o, in terms of ci,3, 03,1 02,2, 02,3, 03,2, and 03,3. (The solutions can be found at the end of the 
section 4.2.) 

Substituting these solutions into our four orthogonality equations gives us four equations and 
six unknowns. These, admittedly long, equations are 

20C3,3C1,3 + 20C3,3C3,1 - 4c2,2Cl,3 - 803,302,2 + 8c3,iCi,3 - 4C3,iC2,2 - 8C3,3C2,3 - 4C3,iC2,3 

-4ci,3C3,2 - 8C3,3C3,2 - 4c3,iC3,2 - 2ci,3 - 2C2,2 - 2c3,i - 2c3,3 + 4c|,2 + 20c|,3 + 8Ci,3 

5 

+4c|2 + 4ci^3 + 8ci,i + - - 4C2,3C1,3 = 4 

-C3,3C3,1 - ci,3 + C3,3 - 03,301,3 - 63,163,2 - 203,363,2 + ^63,2 - 61,363,2 + 6|,2 - 63,162,3 
-263,362,3 + 62,3 + ^62,3 - 62,361,3 - 63,162,2 - 263,362,2 + 62,2 - C2,26l,3 = 

-263,361,3 - 1063,363,1 + 262,261,3 + 463,362,2 - 463,161,3 + 263,162,2 + 463,362,3 + 263,162,3 
+261,363,2 + 463,363,2 + 263,163,2 - 61,3 + 62,2 + 63,1 - 63,3 - 26|,2 - 66|,3 - 26|,2 - 26i,3 - 461,1 

+262,361,3 + ^ = 
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-10C3,3Cl,3 - 2C3,3C3,i + 2c2,2Cl,3 + 4c3,3C2,2 " 4c3,iCi,3 + 2c3,iC2,2 + 4C3,3C2,3 + 2C3,iC2,3 
+2C1,3C3,2 + 4C3,3C3,2 + 2C3,iC3,2 + Ci,3 + C2,2 - C3,i - €3,3 - 20^2 " ^4,3 " ^^j^S " 2ci,2 " 2c|3 

+2C2,3Cl,3 + ^ = 

Notice that all of these equations are quadratic in the six variables, so solutions are basically 
intersection points of hyper-spheres. 

We tried to solve three of the four new equations with Maple, but the output contained six 
possible solution types. Four of these possible solution types led to refinement coefficients with 
imaginary components. The other two gave us real refinement coefficients, but it was clear that we 
did not discover all of the solutions because a solution published in [H] was not included. 

We then found a better way to solve the fourteen equations. The method to solve these four 
equations, giving us the eight types of solutions listed in the next section, is described in detail in 
the Maple file Ssolutions .mws. Basically, the idea is to eliminate equations by combining them 
in such a way as to remove quadratic terms. For instance, subtracting the last two of the four 
equations above yields this much-simpler equation: 

-8C3,3C1,3 + 8C3,3C3,1 + 2ci,3 " 2C3,1 - 4c?_3 + 4ci_i = 0. 

This equation can be solved for ci_3, giving us two solutions. One solution is ci^3 = 03^1, and we 
refer to these as type "A" solutions. The other is ci,3 = — 2c3,3 + | — cs,i, and these are called type 
"B" solutions. There are four versions of each solution, as described in the next section and proved 
in the Maple file. 

4.2 Method of Creating Masks 

In this section, we describe how to create any possible mask of refinement coefficients that satisfy 
our conditions. A mask is a matrix of values for the 16 refinement coefficients. To begin, we must 
recognize that two of the coefficients will become parameters in our final solution, as there are only 
fourteen equations. In our solution, the parameters will be 03^2 and 03,3. 

To create a mask, first choose values of 03,2 and 03,3, along with one of the 8 types of solutions 
listed below. Note that there are requirements (described below) in your choice of 03^2 and 03^3, in 
order to get real-valued solutions. Next, compute four of the refinement coefficients in this order: 
C2,3) C2,2, C3^i, and then 01^3. Finally, the other refinement coefficients can be computed. 

In all of the following solutions, /ii = and /X2 = ■ The following discriminants are 

used: 



Ai = -34 + 3203,3^3 - 20V3 + 3203,2^3 - 48ci^3 + 32c3,3 - 32c3,3C3,2 + 48c3,2 - 48c|2 
A2 = -48ci_2 - 3203,2^3 - 32c3,3C3,2 + 48c3,2 - 34 - 3203,3^3 + 20^3 - 48ci,3 + 32c3,3 

A3 = -48^2 - 32C3,3C3,2 - 48c|3 + 2 - 16C3,3 
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Here are the eight types of solutions: 
Type Ala 



1 1, , 1 ^ 

C2,3 = ^ + 2/Xl--(C3,2 + C3,3) + -VAl 

3 , 

C2,2 = - + 4/Xi - C3,2 - C3,3 - C2,3 
C3,l = fJ-l- C3,3 
Cl,3 = C3,i 



Type Alb 



C2,3 = 7 + 2/xi - 



C2,2 
C3,l 
Cl,3 



+ 4^1 - 

- C3,3 
C3,l 



1/ N 1 

2(C3,2+C3,3)- g- 

C3,2 — C3,3 — C2,3 



Requirement on 03,2 and 03,3 for real solutions of type Ala and Alb: Ai > 0. 
Type A2a 



C2,3 

C2,2 

C3,l 
Cl,3 



- + 2/X2 - -(C3,2 + C3,3) + g V A2 

3 , 

- + 4/^2 - C3,2 - C3,3 - C2,3 
M2 - C3,3 

C3,l 



Type A2b 



C2,3 

C2,2 
C3,l 
Cl,3 



- + 2/.2- 

3 , 
4+4/X2- 

/^2 - C3,3 
C3,l 



1/ N 1 

2(C3,2 + C3,3) - gVA2 
C3,2 — C3,3 — C2,3 



Requirement on 03,2 and 03,3 for real solutions of type A2a and A2b: A2 > 0. 
Type Bla 

C2,3 = -o(C3,2 + C3,3) + ^ V A3 



1 



C2,2 



^ — C3,2 — C3,3 — C2,3 
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Type Bib 



Type Bla 



C3,l = Ml - C3,3 
1 

Cl,3 = 2 ~ ^^^3,3 - C3,l 



1/ ^ nr- 

C2,3 = -2(^3,2 + C3,3) - g V^3 
1 

C2,2 = ^ — C3,2 — C3,3 — C2,3 

C3,l = Ml - C3,3 

Cl,3 = o ~ ^"^3,3 - C3,l 



1/ ^ 1 

C2,3 = -2(^3,2 + C3,3) + g V^3 

1 

C2,2 = ^ — C3,2 — C3,3 — C2,3 

C3,l = M2 - C3,3 
1 

Cl,3 = 2 ~ ^"^3,3 - C3,l 



Type Bla 



C2,3 

C2,2 
C3,l 

Cl,3 



1/ N 1 

-2(C3,2 +C3,3) - gV^3 

1 

^ — C3,2 — C3,3 — C2,3 

M2 - C3,3 

7: - 2C3,3 - C3,l 



Requirement on 03^2 and 03^3 for real solutions of type Bla, Bib, B2a, and B2b: A3 > 0. 

With these six values, the other ten can then be calculated. The following equations are the 
solutions to the 10 linear equations in terms of ci^3, 02,2, £2,3, 03^1, 03^2 and 03^3. 



C0,0 = -C3,l - 2C3,3 - Ci,3 + C2,2 

C0,1 = -C3,l - 2C3,3 + ^ + C2,3 - Cl,3 
1 

C0,2 = Ci,3 + C3,3 + 2 ~ "^2,2 
C0,3 = C3,3 - C2,3 + Ci,3 

Cl,0 = -C3,l - 2C3,3 + C3,2 - Ci,3 + ^ 
11 



Cl,l = -Ci,3 - C3,i - C3,3 + 1 

1 

Cl,2 = -C3,2 + Cl,3 + C3,3 + - 

1 

C2,0 = C3,i + C3,3 - C2,2 + ^ 
1 

C2,l = C3,i + C3,3 + 2 ~ "^2,3 
C3,0 = C3,i + C3,3 - C3,2 



An example of a bivariate Daubeehies scaling function that can be found in the literature is in 
[H], a paper published in 1999 by Bin Han of Princeton. Han's scaling function is of type B2A, 
using C3,2 = H2 and 03,3 = 0. 

In the next section, we will prove an important lemma. Following that is a description of how 
to generate a picture of a scaling function from a mask. Then, we will present masks of all eight 
types, along with the graphs of the related scaling functions. 



5 An Important Lemma 

In this section, we prove another Lemma that will be useful in the rest of this report. 
Lemma 3: ^ • Y,j - hV - j) = 1- 
Proof of Lemma 3: We will first show that 

C2n,2m = X] X] C2n+l,2m = XI X] '^2n+l,2m+l = X X '^2n,2m+l = 1- 

n m n m n m n m 

Define the following values: 

n m 

= Co,0 + Co,2 + C2,0 + C2,2 

'i^2 = X X '^2n+l,2m 
n m 

= Ci,o + Ci,2 + C3,o + C3,2 

= X X/ '^2n+l,2m+l 

n m 

= Cij + Ci,3 + C3,i + C3,3 

Wi = X X '^2n,2m+l 
n m 

= Co,l + Co,3 + C2,l + C2,3 

The first regularity equation can be written as W4 — wi = 0. Similarly the second and third 
regularity equations can be written as wi — = and 1^2 — W4 = 0. Our averaging equation can 
be written as 1^1 + 1112 + 103 + 104 = 4. We now have a system of four equations with four variables. 
The unique solution to this system is wi = W2 = = W4 = 1. 
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Now we will show that Yli Ylj 4>{x — i,y — j) = 1- Define 

f{x,y) = ^^(Pix -i,y - j). (31) 

i 3 

Notice that f{x, y) is continuous, as it is defined in terms of the scaling function which is continuous. 
Using the dialation equation f{x,y) can also be written as: 

f{x,y) = ^^{J2^Cm,n(t>{'^ix-i)-m,2{y-j)-n)) 

i j m n 

= XI XI XI X] Cm,n^(2a; -2i-m,2y - 2j - n) 

i j m n 

Let a = 2i + m and b = 2j + n so m = a — 2i and n = b — 2j. Thus, 

f{x,y) = ^^^^Ca-2i,b-2j(t>{'2x-a,2y-b) 

i j a b 

= X XE X Ca-2j,6-2j)0(2x -a,2y- b) 

a b i j 

We can see that ^a—2i,b—2j — 1 regardless of a and b because of our argument above. Thus, 

a b 

= /(2a;, 2y) 



In order to show that f{x,y) is constant, a proof by contradiction will be used. Assume that 
/(0, 0) = A and for some p and q, f{p,q) = B such that B. Since f{x,y) = f{2x,2y) we can 
deduce that /(|, |) = B, /(|, |) = B, etc. This pattern continues and since f{x,y) is continuous, 
it must be that /(0, 0) = B which is a contradiction of the assumption that A B. Therefore, 
f{x,y) must be constant, let's call this constant P. 

Now we simply need to show that P=l. Since P is a constant, 

1 

Pdxdy 

-^0 

1 /•! 

f{x,y)dxdy 





From the definition of /(x, y). 



Jo Jo i j 

XX / / ^{x-i,y-j)dxdy 
i J Jo Jo 



The previous step is valid because all but a finite number of terms are zero . Using a change of 
variables, we now have 



13 




/ / (l){x,y)dxdy 



J — oo J — oo 



= 1 



by the averaging assumption (0). Therefore, we have shown that (t>{x — i,y — j) = 1. 



Due to the compact support condition for we can actually state that 0(2, 2) + (/>(!, 1) + 
(/.(2,1) + </.(!, 2) = 1. 

6 Generating y) with the Cascade Algorithm 

6.1 Overview 

Once a mask is determined by the analysis explained in the previous section, there is a question of 
how to create the associated scaling function (j){x-, y)- To do so requires a modification of Daubechies' 
cascade algorithm [D]. For the bi-cascade algorithm, we take advantage of the compact support 
condition, which allows us to state that (f){x, y) = when x and y are integers, except at these four 
key points: (1, 1), (1, 2), (2, 1), and (2, 2). We will also make use of the dilation equation (|l|). 

Basically, there are two parts to the bi-cascade algorithm. First, determine the values (p{l,l), 
(/)(1,2), (j){2,l) and (j){2,2) by using a fixed-point method. Then, determine the values at points 
(x, y), where x and y are rational numbers whose denominators are powers of 2. Next, we describe 
each part in detail, and Maple code for the algorithm can be found in the file bicascade .mws. 

6.2 Determination of four key function values 

Due to the compact support condition, the dilation equation (|l|), evaluated at our four key points, 
yields the following four equations: 



Another way to look at this is that if 6 = [(p{l, 1), (p{2, 1), </>(!, 2), (p{2, 2)]^, then we seek solutions 
of the equation b = Lb, where L is the 4-by-4 matrix 



Notice that any vector b which satisfies the equation is an eigenvector corresponding to eigenvalue 
1. 



'^(1,1) 
'^(2,1) 
</>(!, 2) 
<A(2,2) 



ci,i0(l, 1) + co,M2, 1) + ci,o0(l, 2) + 00,00(2, 2) 

C3,10(l, 1) + C2,10(2, 1) + C3,O0(1, 2) + C2,O0(2, 2) 
Ci,30(l, 1) + Co,30(2, 1) + Ci,20(l, 2) + 00,20(2, 2) 
C3,30(l, 1) + C2,30(2, 1) + C3,20(l, 2) + C2,20(2, 2) 



L = 



/ Ci,i Co,l Ci,o Co,0 \ 

C3,l C2,l C3,o C2,o 

Cl,3 Co,3 Ci,2 Co,2 

\ C3,3 C2,3 C3,2 C2,2 / 
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It turns out that one of the eigenvalues of L needs to be 1, while the rest have absolute value 
less than 1 . Recall that the averaging and regularity equations formed a linear system. This system 
determines values for the ten following variables ci^2; C2,i, 02,0, £0,2, co,i, ci^Oi cs^O) ^0,3, and co,o 
in terms of ci^3, 03^1 02,2, £2,3, 03^2, and 03^3. Substituting those values into L and computing the 
eigenvalues with Maple yields four eigenvalues: 1, ^ twice, and —03^2 — £2,3 + C2,2 + £3,3. (See the 
Maple file eigenL.mws.) To show that this last eigenvalue has absolute value less than one, we 
can use Lagrange multipliers, using the orthogonality equations to get two constraints. Subject to 
these constraints, we conclude that 

I - C3,2 - C2,3 + C2,2 + £3,3] < — < 1. 

To determine b, we set up an iterative scheme 

bn+l = Lbn. (32) 

The implication of the analysis of the eigenvalues of L is that this scheme will converege. Further, 
in the original cascade algorithm, all function values of (j) on the integers are initially set to zero, 
except that one of them is set to 1. This way, the sum of the values of <j) on all of the integers will 
be 1. Analogously, for the bi-cascade algorithm, we should initially use 60 = 

[1,0,0,0]'^, and this 

is consistent with Lemma 3. 

Using Maple, the iterative scheme ( 02|) will converge after about 40 iterations. We then have 
the value of (p at the four key points. 

6.3 Determination of other values of (p 

Once the value of (p is known at points where x and y are integers, the values where x and y are 
multiples of ^ can be determined by using the dilation equation (|^). For example, 

i j 

This process can then be iterated to determine the values of (j) on multiples of |, |, etc. As the 
iterations continue, the number of points where (j) is determined increases exponentially, and after 
determining values on multiples of the amount of time needed for the next iterate is not worth 
the effort. 

At the end of the Maple file are several methods for creating graphs of the scaling function. 
One method is a simple plot of points. The second method is to define a function which agrees 
with (j) on all of the determined points. This function is useful for creating linear combinations of 
translates of (j), and will be used in a later section. The third method creates an animated gif that 
rotates the graph of the scaling functions about the z-axis. 

7 Examples of Masks and Associated Scaling Functions 

The following masks were created using the method described in section 4. Each mask is presented 
in the following matrix form: 
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/ Co,0 Co,l Co,2 Co,3 \ 

Cl,0 Cij Ci,2 Ci,3 

C2,0 C2,l C2,2 C2,3 

V C3,o C3,i C3,2 C3,3 / 



In our examples, different parameters were used for the Al examples tlian for the other examples. 
There is no pair (03^2; C3^3) that will yield real valued masks for all 8 types. (See Section 4.2.) 

Ala, using (03,2,03,3) = (1,5): 



-3 - - 46 + 28V3) 

-# + iV3 



|y(-46 + 28V3) 
1 - iv^ 



iV3 + |\/( - 46 + 28^3) 



iv^- i\/(- 46 + 28V3) i + i%/3- i\/( - 46 + 28^3) 
-i+^V3 1 



- jV^- - 46 + 28v^) \ 

-i + iv/3{3) 
|V3+ W(- 46 + 28%/3) 



Alb, using (03,2,03,3) = (l,i): 



-i + iy(- 46 + 28v/3) 
i - iV3- i\/(- 46 + 28V^) 



-ly(- 46 + 28V3) 

[V3+ J7(- 46 + 28V^) 
-i + i^/3 



1 - i%/3- ly(-46 + 28y3) 

i + iV3+ iy(- 46 + 28V^) 
1 



- ly(-46 + 28y3) N 

rv^- 3y(-46 + 28V^) 



A2a, using (03,2,03,3) = (0,0): 



f 3 - 5a/(- 34 + 20v/3) 



I + Iv^C - 34+ 20^3) -i + lv^+ -34 + 20\/3) 

1 + i%/3+ iyi- 34 + 20^3) i\/3- iy(- 34 + 20v^) 1 - i\/3 - iyi - 34 + 20^3) 



H iv^- -34 + 20\/3) ^ 

|a/3+ -34 + 20v/3) 

y 



A2b, using (03,2,03,3) = (0,0): 

i + iyi- 34 + 20 V^) 



- 34 + 20v/3) -i + 1%/3 - - 34 + 20\/3) 

iV3 i + iV3 f-iV3 

i + jV^- IVI - 34 + 20V3) iV3+ - 34 + 20s/3) 1 - i\/3 + W( - 34 + 20a/3) 



^^/3 



^\/3 



-| + iv^+ -34 + 20v/3) ^ 

4 4^3 
|-iV3-iy(- 34 + 20V3) 

/ 



Bla, using (03,2,03,3) = (0,0): 



1V2 







V 



jVs + lV^ 4-i\/2-j\/3\ 



4 V " I 8 
4 8 V ^ 









Bib, using (03,2,03,3) = (0,0): 



V 





5 + 7\/3 



-IV2 
1 
2 



- |\/3 - |\/2 
i + ^V2 



-IV2 

/ 
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B2a, using (03,2,03,3) = (0,0): 



IV2 



\ 



IV3 



IV2 

1 
2 



V2 



+ 



|^/3 + i\/2 



+ 



IV3 




B2b, using (03,2,03,3) = (0,0): 





1 

4 " 



1V2 



iV3 



.1^2 
1 
2 



+ 



+ 



l^-i 4 





8 Reproducing Linear Functions 

8.1 A General Formula to Reproduce Linear Functions 

Just like 1)4, we can reproduce linear functions by using linear combinations of (p{x,y) and its 
translates. To begin, suppose 17 is a function that we wish to reproduce, that is, 

9{x, 2/) = X] X] ^^^'"^i^ -u,y -v) 

U V 

With a change of variables, we can rewrite this equation as 

9{x, y) = X] X] (^x-i,y-j4>{h 3) (33) 

So, using integer values for x and y and an arbitrary bivariate (j) mask, 

g{x, y) = aa._2,y-2'?!'(2, 2) + aa._i,y_i0(l, 1) + a^„2,j/-i0(2, 1) + a:c_i,y_2(/>(l, 2) (34) 

Notice that we only have (/)(2, 2), 0(2, 1), 0(1, 2),and 0(1,1) left because (p{i,j) = at all other 
integer points. 

We are now ready to figure out ax~i,y~j in terms of This will enable us to create linear 

functions out of a bivariate scaling function and its translates. We assume that g{x, y) is a 
polynomial, SO (111 V ^ 

function of u and v will have the same degree as g{x,y). 

If g{x, y) = X, then ai^j = axi + Pxj + 1x for some constants ax, (3x, and 7^.. We will determine 
these constants such that we can write x in terms of 0(2, 2), 0(1, 1), 0(2, 1) and 0(1, 2). Substituting 
into (jsj), and using Lemma 3: 

X = aa;_2,j/-20(2,2) +aa;_i,y_i0(l,l) +a^_2,y-i0(2,l) + a2,_i,y_20(l,2) 

= {ax{x - 2) + (3x(.y - 2) + 7x)0(2, 2) + (q^(x - 1) + ^^(y - 1) + 7^)0(1, 1) 

+{ax{x - 2) + pxiy - 1) + lx)4>{2, 1) + Mx -l)+Px{y-2)+ 7x)0(l, 2) 
= xa,(0(2, 2) + 0(1, 1) + 0(2, 1) + 0(1, 2)) + y/3.(0(2, 2) + 0(1, 1) + 0(2, 1) + 0(1, 2)) 

+a,(-20(2, 2) - 0(1, 1) - 20(2, 1) - 0(1, 2)) + /3,(-20(2, 2) - 0(1, 1) - 0(2, 1) - 20(1, 2)) + 

7.(0(2, 2) + 0(1,1) +0(2,1) +0(1, 2)) 
= xax + yPx + «x(-20(2, 2) - 0(1, 1) - 20(2, 1) - 0(1, 2)) 

+/?,(-20(2, 2) - 0(1, 1) - 0(2, 1) - 20(1, 2)) + -fx 
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Therefore Ux = I and Px = 0, and 7^ = 20(2, 2) + (f){l, 1) + 2(p{2, 1) + 0(1, 2). Using these constants 
we get 

X = ^ ^(i + 20(2, 2) + 0(1, 1) + 20(2, 1) + 0(1, 2))0(x -z,y~ j) (35) 
i i 

Similarly, we will write y in terms of 0(2, 2), 0(1, 1), 0(2, 1) and 0(1,2). If g{x,y) = y, then 
aj j = Oyi + Pyj + 7y for some constants Oy, (3y, and jy. Again using ( p^ and Lemma 3 gives us: 

y = aa-2,j/-20(2, 2) + aa._i,y_i0(l, 1) + aj;_2,j/-i0(2, 1) + 0^.-1,^-20(1, 2) 

= {ay{x - 2) + Py{y - 2) + 7^)0(2, 2) + {ay{x - 1) + /3y(y - 1) + 7^)0(1, 1) 

+{ay{x - 2) + /3y(y - 1) + jy)^{2, 1) + (aj,(x - 1) + Py{y - 2) + 7j,)0(l, 2) 
= xaj,(0(2, 2) + 0(1, 1) + 0(2, 1) + 0(1, 2)) + 2//3j,(0(2, 2) + 0(1, 1) + 0(2, 1) + 0(1, 2)) 

+a^(-20(2, 2) - 0(1, 1) - 20(2, 1) - 0(1, 2)) + /3j,(-20(2, 2) - 0(1, 1) - 0(2, 1) - 20(1, 2)) + 

7^,(0(2, 2) +0(1,1) +0(2,1) + 0(1, 2)) 
= xay + ypy + aj,(-20 - 0(1, 1) - 20(2, 1) - 0(1, 2)) + 

Py{-2(t> - 0(1, 1) - 20(2, 1) - 0(1, 2)) + Jy 

Therefore ay = 0, Py = 1, and jy = 20(2, 2) + 0(1, 1) + 0(2, 1) + 20(1, 2). Thus, we have found 
constants ay,Py, and 7^ such that g{x,y) = y, and 

y = Yl + 2) + 0(1, 1) + 0(2, 1) + 20(1, 2))0(:e -i,y- j) (36) 

i j 

In equations (^) and (|36[), along with Lemma 3, we now have x, y and 1 written in terms 
of 0(2, 2), 0(1, 1), 0(2, 1), and 0(1,2) . This enables us to write an approximation for any type of 
plane given a particular mask. 

Let us look at the general case of f{x, y) = kx + ly + m where k, I, and m are constants. Thus, 

f{x,y) = kx + ly + m 

= KY. + 2) + 0(1, 1) + 20(2, 1) + 0(1, 2))0(x - i, y - 3)) 

i j 

+ 2) + 0(1, 1) + 0(2, 1) + 20(1, 2))0(x -i,y- j)) 

i j 

- hy - j)) 

i 3 

= + 20(2, 2) + 0(1, 1) + 20(2, 1) + 0(1, 2)) + l{j + 20(2, 2) + 0(1, 1) + 

« j 

0(2,l) + 20(l,2))+m)0(x-i,y-j) 
This simplifies to 

kx+ly+m = Y ^(2-^(2, 2)(/c+/)+0(l, l)(A:+/)+0(2, 1)(2A;+O+0(1, 2){k+2l)+ki+lj+m)(l){x-i, y-j) 

i j 

(37) 

Using the above equation, we can generate f{x, y) using a bivariate scaling function associated with 
a given mask. 
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8.2 An Example of Reproducing A Lineeir Function 

Consider the function g{x, y) = x + 6y — 10, a linear function that we wish to reproduce using one 
of our scaUng functions and its translates. In this case, k = 1,1 = 6, and m = —10. So, 

9i^^y) = EE(2(1) + 2(6M2,2) + ((1) + (6M1,1) + (2(1) + (6M2,1) 

i 3 

+((1) + 2(6))<^(1, 2) + {l)i + (6)i + (-10))<^(x -i,y- j) 
= E 2) + 7(/.(l, 1) + 80(2, 1) + 130(1, 2)+i + 6j- lO)0(x -i,y- j) 

i j 

We will reproduce g using the scaling function in the figure below. This is a B2bl type of scaling 
function, when 03,2 = £3,3 = 0. Note that this scaling function is also plotted in section 4.2. 

/ -I-IV2 IV2 1 + 1^3 + 1^/2 1-IV2 + IVS\ 

i l + \Vs i + iV3 

i- 1^/3 + 1^/2 j-jV^-^V^ iV2 

\ l-lVs \-lVs / 

9 Areas for Further Study 

One area for further study is exploring the mother wavelets. (Once a scaling function is defined, 
a mother wavelet can also be defined. See [AS] and [D].) According to Daubechies, there should 
be at least three for each of our bivariate scaling function 0. Another idea is to develop a further 
classification scheme for these scaling functions. A bivariate version of Dq would also be interesting, 
as they could be used to exactly reproduce bowl and saddle surfaces. Also, our undocumented 
investigations into bivariate Haar wavelets can be tied into the work in this report. 
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11 Note about JVIaple files and Image Files 

A longer version of this paper is available which includes image files of surface graphs of bivariate 
scaling functions. It can be downloaded at: 

www . gvsu . edu/mathstat/wavelets/undergrad . htm 

Maple files may be requested from Edward Aboufadel. 
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